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ABSTRACT 

In this paper we present emission line strengths, abundances, and element ratios 
(X/0 for Ne, S, CI, and Ar) for a sample of 38 Galactic disk planetary nebulae (PNe) 
consisting primarily of Peimbert classification Type I. Spectrophotometry for these PNe 
incorporates an extended optical/near-IR range of AA3600-9600 A including the [S III] 
lines at 9069A and 9532A, setting this relatively large sample apart from typical spectral 
coverage. We've utilized ELS A (Emission Line Spectrum Analyzer), a 5-level atom 
abundance routine, to determine Te,Nc, ICFs, and total element abundances, thereby 
continuing our work toward a uniformly processed set of data. With a compilation 
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of data from >120 Milky Way planetary nebulae (PNe), we present results from our 
most recent analysis of abundance patterns in Galactic disk PNe. With a wide range of 
metallicities, galactocentric distances, and both Type I and non-Type I objects, we've 
examined the alpha elements against HII regions and blue compact galaxies (H2BCG) 
to discern signatures of depletion or enhancement in PNe progenitor stars, particularly 
the destruction or production of O and Ne. We present evidence that many PNe have 
higher Ne/0 and lower Ar/Ne ratios compared to H2BCGs within the range of 8.5-9.0 
for 12-|-log(0/H). This suggests that Ne is being synthesized in the low-intermediate 
mass progenitors. Sulfur abundances in PNe continue to show great scatter and are 
systematically lower than those found in H2BCG at a given metallicity. Although we 
find that PNe do show some distinction in alpha elements when compared to H2BCG, 
within the Peimbert classification types studied PNe do not show significant differences 
in alpha elements amongst themselves, at least to an extent that would distinguish in 
situ nucleosynthesis from the observed dispersion in abundance ratios. 

Subject headings: ISM: abundances, nuclear reactions, nucleosynthesis, abundances, 
planetary nebulae: general, stars: evolution 



Introduction 



Planetary Nebulae (PN) abundance patterns have long been used to note signatures of nuclear 
processing and to trace the distribution of metals throughout galaxies. Being the shed envelopes 
resulting from the late evolutionary stages of intermediate-mass stars (IM© < M < SMq), PNe 
abundances reflect both the evolution of their progenitors and the natal interstellar medium from 
which they formed. PNe are populous, cover a range in progenitor mass and metallicity, and are 
distributed throughout the disk, bulge, and halo of our Galaxy; hence their abundances are useful 
probes into the evolution of intermediate mass stars, the yields of previous generations of stars, 
and the chemical evolution of galaxies. 

This work continues our compilation of abundances in Galactic PNe that are based upon newly 
acquired CCD spectrophotometry over an extended range of optical/near IR wavelengths including 
the strong [S III] emission features at A9069 and A9532 A. Here we present new spectrophotometric 
data for 38 mostly Type I Galactic PNe. These new spectra, wh en added to our previous samples 



(jMihngo et al. I l2002bl : iKwitter et al. I l2003l : iHenrv et al. 1 12004 ) bring the total sample to >120 



covering a substantial range in galactocentric distance and including Peimbert classification Types 
I and II as well as several anticei iter PNe. Peimb ert "types" are categorized by the degree to which 
N and He have been enhanced ([Peimbert Ill978l ). This is indicative of progenitor mass insofar as 
there is a mass dependence on the sites of nucleosynthesis that create and alter the abundances 
we measure in PNe. Type I PNe have high N and He which is consistent with th em having 



progenitors > 2.5 M0, whereas Type II PNe are inferred to have less massive progenitors (jPeimbert 
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19781 : iMaciel I Il992l : iKingsburgh &: Barlow Ill994l ). With this data in hand we are compelled to 
look for distinguishing signatures of nuclear self-contamination across the range of metallicity, 
Peimbert "type", and galactocentric distance our sample contains. The questions of oxygen and 
neon depletion and/or enhancement as seen in the PN gas are of particular interest. 

The evolution of PNe progenitors gives rise to convective episodes that reach down into tem- 
perature regions where nuclear processing can occur. The subsequent dredge-up of products from 
this in situ nucleosynthesis alters the natal envelope material. The first dredge-up (FDU) occurs on 
the first giant branch and is a significant mixing event for the lower-mass end of PNe progenitors < 



1.5 M0 (jBoothrovd &: Sackmann Ill999l ). While on the asymptotic giant branch (AGE) the third 
dredge-up (TDU) and hot bottom burning (HBB) may occur, depending on the mass and metal- 
licity of the progenitor. The ejected envelopes of these evolved stars become planetary nebulae, so 



(van den Hoek & Groenewegen 


1997; 


Renzini &: Voli 


1981 



that Type I PNe have long been connected to bipolar morphology (jPeimbert fc Torres-Peimbert 



It is important to note here 



1983) and that bipolar morphology in PNe is most likely the result of binary progenitors (Ide Marco 



20o3). The effects of binary interactions on the evolution of the progenitor and the composition of 
the subsequent PNe gas could be at play here but are beyond the scope of this paper. 

Within the published literature there is a lack of consistent observational support for the 
production/depletion of oxygen and neon in PN abundances. Oxygen has a significant history 
here as it is a common metallicity tracer for the interstellar medium; marking the degree to which 
chemical enrichment has occurred. However oxygen is problematic in PNe as it is suspected of self- 
contamination, in other words undergoing nuclear processing during the evolu tion of the prog enitor 



star. Oxygen is subject to depletion through the ON cycle during HBB (jClavton I Il968l ). and 
possibly enrichment via the TDU of ^^O from the pr oducts of He burning and subsequent alpha 



capture (jiben &: Renzini Ill983l ; iPequignot et al. Il200d ) . The degree to which the surface abundance 
of oxygen is expected to be altered is unclear as it is subject t o the parameter space of dredge- 
up episodes and HBB (IRenzini k, Voli Ill98ll ; ICharbonnel I l2005l ) . In PNe progenitors ~ 3-5 M0, 
oxygen depletion owing to HBB could be modestly visible though difficult to discern within the 
typical uncertainties associated with PN abundance work (0.1-0.3 dex). Alternate tracers such 
as sulfur and argon, which are believed to be precluded from self-contamination, could instead 
be used to indicate the metallicity of the natal ISM of PN progenitors. Neon also has a history 
as a metallicity tracer but it is also suspected of undergoing self-contamination. Neon can be 
enriched through the production of ^^Ne in TP-AGB stars via two a captures onto ^^N Jbusso 



20061 ; iKarakas &: Lattanzio II2003I ; iMarigo et al. II2003I ) . According to iKarakas h Lattanzio I (|2003l ) 



the production and dredge-up of Ne is particularly efficient in progenitor stars ~ 2-4 M0 across 



( Henrv 


1989: 


Henrv et al. 


2004) 



If we do 

expect elemental abundances to be altered due to self-contamination and subsequent dredge-up, 
then observing patterns that indicate depletion or enhancement across the parameter space of 
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PNe serves as observational support of our understanding of the nuclear processing and convective 
episodes that occur during the late stages of low and intermediate-mass stellar evolution. 

Observationally it has proven difficult to note modest depletions and /or enhancements in alpha 
elements in PN. The scatter seen in observed PN abundances is due to many factors including 
uncertainties in total element abundances, possible systematic effects brought on by compiling 
inhomogeneous data sets, as well as the true spread in elemental abundances that is expected due 
to the range of PN progenitor masses, metallicities, galactocentric distance (abundance gradient 
effects), and the degree of nuclear self-contamination. We are attempting to address some of this 
using extended spectral coverage, allowing for good quality sulfur abundances, and a uniformly 
processed set of data, in other words internally consistent observation, reduction, and abundance 
determination techniques. 

Subsequent sections of this paper discuss data acquisition and reduction, our abundance deter- 
mination scheme, as well as results and analysis of X/H and X/0 ratios for O, Ne, S, CI, Ar, and 
N. An important feature of our work is the use of the near-IR [S III] lines and, where available, a 
[S III] temperature to calculate S"*"^ abundances. Because a significant portion of S in PNe resides 
in and higher ionization stages, including the near-IR [S III] lines improves the extrapolation 
from observed ion abundances to total element abundance. In § 4.2 we discuss issues related to 
observed sulfur abundances in PNe. 



2. Observations and Reductions 

We've acquired spectra for 38 Galactic disk PNe with coverage that extends from 3600A 
to 9600A. Using the same instrumentation and setup as our previous Type II PNe program 



(jMilingo et al. Il2002al ). observations were obtained over a total of four runs; Kitt Peak National 
Observatory (KPNO) in June 2003 & August 2004, and Cerro Tololo Interamerican Observatory 
(CTIO) in November 2003 & August 2004. Table 1 lists the objects, total exposure times in the 
blue and red grating configurations, and the run during which the spectra were acquired. 

To obtain the needed range of spectral coverage, at both observatories we used a combination 
of two gratings with overlap in the region of Ha for subsequent merging. The KPNO observations 
utilized the 2.1m telescope and GoldCam spectrograph. We used a long E-W slit with projected 
dimensions of 5"x 305"to accommodate angularly large objects and simultaneous "sky" exposure. 
Spatial scale for the GoldCam is 0.78 "/pixel. On the GoldCam we used gratings 240 & 58 in first 
order with nominal wavelength dispersions of 1.5 A/pixel and 1.9 A/pixel respectively. This gave 
us a FWHM resolution of 8A and lOA respectively. 

The CTIO observations were made with the 1.5m telescope and the R-C (Cassegrain) spec- 
trograph using a long E-W oriented slit with projected dimensions of 5"x 321". The spatial scale 
for the R-C spectrograph is 1.3 "/pixel. Gratings 9 & 22 were used in first order with nominal 
wavelength dispersion of 2.8 A/pixel and 8.6A FWHM resolution for both. 
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The angular sizes within our Type I sample vary greatly, so we either centered the PNe in the 
slit or concentrated on a bright portion of the nebula. Any offsets from the center of the PNe are 
listed in Table 1. If possible the central star was avoided during the observations and extraction 
process, but for angularly small objects it is unavoidably included. 

Standard calibration frames were obtained for bias subtraction, flat-fielding, slit illumination 
correction, and wavelength calibration. Standard star spectra were also obtained for flux calibra- 
tion. The removal of instrumental signature, calibration, and extraction of the final one-dimensional 
spectra were done using the twodspec package in IRAF[^. At both KPNO and CTIO, the thinned 
red chips produce interference fringes visible in the red. In our spectra the fringes appear at the 
~ibl% level at ~7500A and increase in amplitude with increasing wavelength: ±1-2% at 8OOOA, 
±4-6% at 8500A, ±6-8% at 9000A. Even at their worst, i.e., at ~A9500, the longest wavelength we 
measure, the fringe amplitude reaches only about ±7-10%, and we note this additional uncertainty 
in our line intensities longward of ~7500A. 



3. Line Strengths and Abundances 

Raw line fiuxes were measured in IRAF using the splot task. For each object the red and blue 
ends of the spectrum were merged using the overlap in the Ha region. Table 2 lists the raw and 
dereddened line strengths for our sample of PNe. Columns labeled F(X ) give raw uncorrected fiuxes 
as measured in splot. These F(\ ) values have been scaled to H/3 = 100 using the observed value of 
FfHP) for each object. Intensities of strong lines have measurement uncertainties of <10%, single 
colons indicate uncertainties of >25%, and double colons indicate uncertainties >50%. The last 
three rows of the table give the extinction measure, c, for each object, the final converged values for 
Ha/H/3, and the log of the observed H/S flux in units of ergs cm~^ s~^. Final dereddened intensity 
ratios are calculated as follows: 

i(AL ^ mxiocf(A) 

l(H/3) F(H/3) 

where I(A)/I(H/3) is the final intensity ratio of an emission line to H/3, FA/F(II/3) is the observed 
fiux ratio of the line to H/3, c is the extinction measure, and f(A) is the reddening function. 

To calculate ionic and total element abundances we have e mployed ELSA (Emission Line Spec- 



trum Analyzer), described in detail by iJohnson et al. I (j2006l ). ELS A incorporates a 5-level atom 



abundance program, spectrum merging and de-reddening functions (jSavage &: Mathis Ill979l ). and 
routines for determining values for the interstellar reddening measure (c) and nebular temperature 
and density (Tc and Ne). If there is a non-zero fiux present for He II A4686, ELSA corrects for 
contamination of the Balmer lines due to ionized helium. This iterative correction exploits the 



^IRAF is distributed by the National Optical Astronomy Observatories (NOAO), which are operated by the 
Association of Universities for Research in Astronomy (AURA), Inc., under cooperative agreement with the National 
Science Foundation. 



- 6 - 



temperature and density dependence of the recombination coefficients for hydrogen and hehum, 
ultimately converging on a value for the intrinsic Hct/H^S ratio. ELS A also allows for error propa- 
gation from measured line strength uncertainty; for example, the uncertainty in a final line intensity 
includes the estimated uncertainty in the line flux as well as the uncertainty in c. Uncertainties 
in temperature and density diagnostics incorporate the errors in their constituent lines, and ionic 
abundance uncertainties include errors in the relevant line intensities, temperature, and density. 
Final abundance uncertainties incorporate errors in the ionic abundances. ELSA does not consider 
uncertainty in the ionization correction factors (ICF), nor does it include any systematic sources 
of uncertainty (e.g., flux calibration errors). 

Table 3 lists the electron temperatures and densities calculated for each object. The temper- 
atures have been determined from diagnostic emission line ratios of [O III], [N II], [O II], [S II], 
and [S III]. Electron densities were determined using emission line ratios of [S II] and [CI III] (see 
the footnotes to Table 3 for more detail). Uncertainties in Tg and Ne are propagated from the 
uncertainties in the line fluxes used to calculate them. 

In calculating ionic abundances, ELSA uses the temperature appropriate to the ionization 
stage whenever possible, leading to a two-zone model. The [O III] temperature is used to calculate 
ionic abundances for He"*", He"*"^, O"*"^, Ar"*"^ and higher, Cl"*"^ and higher, Ne"*"^ and higher, and S"^^ 
when T[g mj is not available or differs by more than 5000K from T[o iii]- The [N II] temperature 
is used for all other ions: O'', O"'", N"*", S"*", and CI'''. Table 4 gives the temperature used in 
determining the ionic abundance for each observed transition. 

The general scheme for determining total element abundances is to sum the observed ionic 
contributions and correct for any unseen ionization stages with an ionization correction factor 
(ICF). We symbolically represent this scheme as: 

n(x)-{E^}-'cf(x). 



Here 1^ is the dereddened H/3-scaled value of the line strengths of observed ionic species, the 
factor is the energy generated per transition for a parti cular ion, ICF(X) is the io nization correction 
factor determined using the functional forms given in iKwitter &: Henry I (|200ll ). and finally N(X) 
is the total number abundance of element X with respect to H. 

Table 4 shows the observed ionic species and associated transition wavelengths (in A) in column 
1, their respective number abundances (as determined by that particular transition) in column 3, 
and the temperatures used to determine those abundances in column 2. [S II] densities were used for 
all ionic abundance calculations. Each final ionic abundance is given as the mean of the observed 
transitions (those marked with an asterisk), weighted by flux. The ionization correction factors 
used in calculating the total element ab undances are also shown. As stated before ELSA does not 
consider uncertainty in the ICF. From iKwitter &: Henrv I (|200ll ) the uncertainty in the ICF for 
oxygen is approximately 10% and the ICFs for other elements are uncertain by 20-30%. Table 5 
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shows the total element abundances X/H and element ratios X/0 for each object. 



4. Analysis and Discussion 
4.1. Alpha Elements 



The primary motivation for this project is to test for evidence of alpha element destruction 
or synthesis during the evolution of PN progenitors. The abundances of Ne, O, S, c|§, and Ar 
have been measured in nearly all of the sample objects. Thus it is possible to study their behavior 
relative to one another in PNe, as well as to compare abundance patterns in PNe with what is 
observed in a combined sample of H II regions and blue compact galaxies (BCGs). The latter 
two object types provide probes of relatively homogeneous interstellar material which generally 
show little evidence of abundance enhancements produced by the evolved stars harbored within 
them (jWofford II2009I ). If these alpha elements are produced in lockstep from a given demographic 
distribution of stars, we should see their co-evolution clearly displayed in PNe and H2BCG. If PNe 
progenitors are altering these alpha elements in the late stages of evolution it can potentially show 
up as a departure from the lockstep behavior. 

Figures [llfTOl show plots of the ten element pairs, where the quantity on each axis is the 
logarithm of the number abundance of the ele nient with respect to H, assuming that logII=12. 
We include three Galactic PN samples, that of iHenrv et al. I (|2004l ) of mostly Ty pe II PNe, the 



Type I objects of the c urrent paper, and PNe located in the anticenter direction from lKwitter et al. 



(j2009. in preparation! ). These three samples are referred to in the legends as HKB04 , MilingoOQ, and 



Bernard-Salas et al. 



KwitterOQ, respectively. In addition we have added LMC objects recently studied by[ 
([20081), referred to in the legends as BS, for the discussion of sulfur abundances in PNe. In all, 
the complete PN sample covers a metallicity ra nge of roughly 7.6 t o 8.9 in 12+log(0/H). For 



comparison, we note that the sample analyzed by IWang &: Liu I (120081 ) extend down to 7.3. 



The H II baseline in the figures comprises data of blu e co mpact galaxies fromllzqtov &: Thuan 



(jl999l ) , Galactic H II regions from IShaver et al. I (jl983l ) and lEsteban et al. I (119981). H II region 



data f or external galaxies from lSmith I (|l975l ). and MlOl data from lKennicutt. Bresolin. &: Garnett 
' Collectively these will be referred to henceforth as H2BCG. Finally, abundance data for the 



sun (jAsplund et al. mod ) appear in our plots as well. Symbols are defined in the figure legends. We 
note that CI abundances were not available in the H2BCG data, and thus Figs. El [H O and 1101 show 
PN abundances only. It should be noted that we are making an a priori distinction between Type 



^While CI is technically not an alpha element, since its stable isotopes Cl'^^ and Cl^^ have odd atomic masses, we 
nevertheless include it in our study of alpha elements be cause it arises a s a secondary product during oxygen burning 
in stars when alpha particles are in high concentration ( Clayton 20031 ) . 

^We note that there is a paucity of neon abundance data in high metallicity H II regions. This stems partially 
from the fact that in metal-rich systems, the [Ne III] emissions are often weak, resulting in poor signal-to-noise ratios. 
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I and non-Type I PNe. As discussed in iHenrv et al. I (|2004l ) Type I PNe have enhanced nitrogen 
(and often helium) abundances and are inferred to come from more massive progenitors hence from 
a more recent epoch. Type II PNe have less N and He enhancement and are inferred to come from 
a more intermediate age population with lower progenitor masses than Type I PNe. The original 
Type I classifica tion included both enh a nced He/H and N/0. In this work we continue to follow 
the discussion of iKingsburgh fc Barlow I (jl994l ) who use only nitrogen as a discriminant, setting the 
minimum nitrogen abundance for Type I PNe as the sum of the carbon and nitr o gen a bundances 
seen in the H II regions in the host galaxy. Given this IKingsburgh fc Barlow I (119941 ) calculate 
the minimum N/0 ratio for Type I PNe in the Galaxy to be 0.8. Using mo re recent abundance 



measurements for the Orio n Nebula (lEsteban et al. Ill998l ) and for the Sun (jGrevesse fc Sauval 



1998; 



Asplund et al. II2005I ). iHenrv et al. I (j2004l ) derive a lower minimum N/0 ratio for Type I 



classification. For this paper all PNe with N/0 > 0.65 are classified as Type I. We first discuss the 
establishment of the H II baseline used in the figures before proceeding with our consideration of 
the PN data. 



4-. 1.1. H II Regions and Blue Compact Galaxies: A baseline for comparison 



In Figs. [TIfTOl the most notable feature is the relatively tight linear behavior (in log-log space) 
demonstrated by the II2BCG objects (except in the cases of CI, as noted). The parameters for the 
least squares fits to the H2BCG data, shown with solid lines in the figures, are provided in the top 
section of Table [U where the left column identifies the relation, while columns 2-5 give values of 
the y intercept, slope, correlation coefficient, and number of sample objects included, respectively. 
Quantities in parentheses indicate the number of outliers removed before calculating the fit. Note 
that uncertainties for individual sample objects have not been taken into account. The first six 
relations (the last two will be discussed later) possess very high correlation coefficients, indicating 
the correlations are unlikely to be statistical accidents. 

Results for H2BCGs in Table [6] show that all slopes are consistent with a value of unity, in 
agreement with the conventional wisdom regarding alpha element p roduction. This last point is of 
great interest, since over a slightly larger metallicity range (7.3-8.8) IWang &: Liu I (j2008l ) inferred a 



slope for Ne vs. O of -|-1.14ib.01 for a sample of extragalactic H II regions and concluded that the 
rate of cosmic Ne evolution increases with metallicity. However, we find that the average log(Ne/0) 
for objects with 12-|-log(0/H)<8.5 is -0.70 while the average for objects above this threshold is -0.75, 
i.e., within the uncertainties the two sets are indistinguishable regarding Ne/0. 



Finally, we employed the chemical evolution code developed by R.B.C. Henry [see lHenrv &: Prochaska 



(j2007l ) for a detailed discussion of the code] to predict interstellar abundances of O, Ne, S, CI, and Ar 
as a function of time over the age of the Galaxy. We then formed theoretical element to element plots 
like those in Figs.fll fTOl calculated slopes, and compared the slopes wit h their observed counterparts 



The m odels included the massive s tar yields for O, Ne, and S from iPortinari. Chiosi. &: Bressan 



(| 19981 ;) and CI and Ar yields from Iwooslev fc Weaver I (Il995l l and SNIa yields of 



Nomoto et al 
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(|l997l ). while the yields of low and intermediate mass stars taken from iMarigo I (|200ll ) do not 
include these elements. Thus, we are only testing the influences of massive stars and SNIa on 
alpha element behavior. Table [7] provides a comparison of the model-predicted slopes with those in 
column 3 of Table [H where we see good consistency between observation and theory, especially in 
light of the large uncertainties usually associated with stellar yield predictions. However, because 
of the yield uncertainties, we are unable t o use the models to t est whether there is a small increase 
of Ne/0 with metalhcity, as suggested by IWang fc Liu I (120081 ). 



We note that the sun's position in Fi g. [T] falls below the H2 BCG trend line by about 0.15-0.20 
dex, according to t he solar abundances of I Asplund et al. I (|2005l ). This offset is consistent with the 
recommendation of I Wang &: Liu I (|2008l ) to raise the solar Ne value by 0.22 dex. Since we conclude 
in the next subsection that Ne is enhanced in PNe, we would not recommend raising the solar Ne 
value to an extent which would place it above the trend line. 

In summary, the evidence found in H2BCG objects appears to be consistent wit h the expected 



lockstep behavior among the alpha elements. This conclusion is at odds with that of I Wang &: Liu 
(|2008l ) in the case of Ne versus O. 



4-1-2. Planetary Nebulae 

We now proceed to study PN abundance patterns with respect to the established H2BCG 
baselines contained in Figs. [TlllOl In doing so, we note that three of the four plots involving sulfur. 
Figs. [21 [5l and [U display what we have previously described in IIKB04 as the sulfur anomaly, 
i.e., the systematic tendency for PN sulfur abundances to be less than those of II2BCG at the 
same metalhcity. This topic is addressed in a separate subsection. We presently focus on the PN 
distribution in the remaining six plots and make comparisons with the abundance distributions of 
the II2BCG objects where possible. 

We carried out a linear regression analysis on the PN data in Figs. [TlllOl for the six cases not 
involving sulfur. The resulting parameters are provided in the bottom portion of Table [6l While 
in the first three cases, Ne-0, Cl-0, and Ar-0, the slopes are consistent with unity, the slopes for 
Cl-Ne, Ar-Ne, and Cl-Ar are not. In fact the slopes of the CI vs. Ne and Ar vs. Ne relations 
differ significantly from unity. If we compare the slopes of PNe and H2BCG in Table 6 for the 
same relations, we see that agreement is present in the cases of Ne-0 and Ar-0 but not in the 
case of Ar-Ne. Here the PN slope (-1-0.57) is much shallower than the slope for H2BCG (-1-0.96). 
We suspect that the flatter slope for PNe is due to the six low metalhcity PNe residing above the 
H2BCG track in Fig. [3 By excluding these objects one sees that the slope would be closer to 
unity. The slopes of the other two element relations which differ from unity, Cl-Ne and Cl-Ar, are 
likely influenced by the greater uncertainty in the CI abundances, given the weak line strengths 
of the CI ions in our spectra. Thus, after this initial look PN abundance patterns appear to be 
reasonably consistent with the expectations of lockstep behavior predicted by theory. If any of the 
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alpha elements are enhanced or depleted by PN progenitor stars, the effect is too small to manifest 
itself in these plots. We now take a closer look at O, Ne, and Ar. 

O, Ne, and Ar: We searched for potential differences between PNe and H2BCG by forming 
ratios of Ne/0, Ar/Ne, and Ar/0, and then comparing the PN and H2BCG frequency distributions 
in each of the three instances. Such histograms could potentially reveal significant differences 
between the two object types if nuclear processing is occurring in the PN progenitors. The resulting 
histograms are displayed in Figs. [TT]|131 The numbers in the legend give the average and standard 
deviation of the log of the relevant element ratio. 

Fig. [TT] shows that Ne/0 in PNe is slightly higher on average than in H2BCGs although the 
difference in the averages is only 0.06 dex, much smaller than the dispersion in the values for PNe. 
More importantly, while the H2BCG distribution is symmetrical, the PN distribution is not, being 
noticeably more extended in the direction of higher values of Ne/0 and consistent with the tendency 
in Fig. [1] for PNe to fall above the H2 BCG trend line. Als o shown in Fig. [11] are distributions for 



the LMC and SMC PNe reported in iLeisv fc Dennefeld I (|2006l ). Since the order of increasing 



metallicity for the four samples is SMC, LMC, H2BCG, MW disk PNe, there is an indication her e 



that Ne/0 in PNe increases with metallicity, in agreement with the findings of lWang &: Liu I (|2008l ). 
While temperature-sensitive lines such as [O III] A4363 generally become weaker and more difficult 
to measure as metallicity rises (due to more efficient cooling) a simple numerical exercise reveals 
that the observed positive offset in Ne/0 from the H2BCG line would require that temperature be 
in error by several thousand degrees below the real temperature. Since in general weak lines are 
usually measured to be stronger than they really are, that would push A4363 and the corresponding 
[O III] temperature up, i.e., in the opposite direction necessary to explain the high values of Ne/0. 
We can therefore eliminate temperature measuring errors as the source of the systematically high 
Ne/0 derived for PNe. 

The average value for log(Ar/Ne) in PNe, according to Fig. [T2t is less than that of H2BCGs by 
0.14 dex, and again the distribution is broader for PNe, with the extension this time appearing on 
the low side. In this case the offset may be significant, since it is roughly one standard deviation in 
magnitude. Finally, Fig. [13] indicates that PNe on average have slightly lower values of log(Ar/0), 
although the difference is very small. At the same time, we note that now the PN distribution 
appears symmetrical. 

Thus, compared with H2BCGs, Ne/0 in PNe seems to extend to higher values, Ar/Ne extend to 
lower values, and the Ar/0 distribution appears very similar. Taken together these three histograms 
suggest that some Ne enhancement is present in a significant portion of PNe, while there is no 
evidence for similar shifts in Ar and O. If correct, this is an important result, as we can now say 
that the elevated Ne/0 ratio in PNe is due specifically to Ne enhancement and not O depletion. 
Furthermore, it suggests that LIMS, the stellar progenitors of PNe, actually produce some Ne 
during their evolution. 



One of the primary motivations for this study was to look for differences between Peimbert 
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Type I and non-Type I PNe in terms of alpha element nucleosynthesis in their progenitor stars. To 
explore this question, we considered the offsets of PNe from H2BCG samples for the two element 
pairs Ne-0 and Ar-Ne. We first computed Ne/0 and Ar/Ne for all objects, determined a least 
squares fit to the H2BCG ratios in both cases (reported in the last two entries of the top panel in 
Table[6|) and then calculated PN offsets, i.e., d(Ne/0) and d(Ar/Ne). In each case, for element ratio 
X, d(X) = XpN — Xh2BCG- We then separated the PN sample into Type I and non-Type-I objects 
to see if there is a difference between these types in terms of the offsets d(Ne/0) and d(Ar/Ne) from 
H2BCGs. Figs. [H] and [E] display these comparisons. In both figures the two PN types appear 
to overlap closely and their averages differ by less than a standard deviation. Thus, we see no 
compelling evidence to suggest that Type I and Type II progenitors differ with respect to their 
distributions o f their offsets fro r n the II2BCG sample. This conclusion is supported theoretically by 
recent work of iKarakas et al. I (j2008l ). who find that Ne is only produced in low and intermediate 
mass stars over a very narrow mass range. A significant difference between PN types is not expected 
in regards to alpha element nucleosynthesis. 

In conclusion, we have presented evidence to support the idea that many PN progenitors 
possess enhanced levels of Ne compared to H2BCG, although O and Ar appear to remain unchanged 
during the star's evoluti on. If this Ne enhancement is real it is ma rking in situ production of ^^Ne 



and the TDU in action (Karakas fc Lattanzio 



2003 



Busso II2OO6I ). However, we see no difference 



between Type I and non-Type I PNe in terms of alpha element nucleosynthesis, despite the fact 
that the former exhibit enhanced N/0 and He/H which are familiar empirical signatures of nuclear 
processing in LIMS. This point is emphasized by our finding that Type I and Type II PNe, which 
differ greatly in N/0, demonstrate no apparent difference in Ne/0. 

We note the presence of a small group of outliers especially in the case of Ne and O in Fig. [TJ 
Their deviation from the norm indicates that in some cases either alpha destruction or production, 
particularly in the case of Ne, does in fact occur, or alternatively that these objects were formed 
out of unmixed interstellar materials such that their composition does not accurately sample their 
current mixed environments. Furthermore, we caution that the presence of those same outliers may 
help to hide systematic trends such as an increase in Ne with metallicity. After all, it is difficult 
to decide the guidelines for labeling an object an outlier and subsequently removing it from the 
sample when we are searching for systematic trends. It is possible that subtle trends connected with 
alpha element evolution do exist in PNe but are concealed by the uncertainties in our abundance 
determinations or the presence of outliers. 

Finally, we find that Ne/0 is systematically higher in the general PN population compared 
to H2BCGs within the range of 8.5-9.0 for 12+log(0/H). This is consistent with the finding of 
Wang fc Liu1 (|2008l v 
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4.2. The Sulfur Anomaly 

We noted previously that PN sulfur abundances appear to be significantly offset from the 
H2BCG track in Figs. [21 [5l and[9j In addition, the scatter in S abundances for a given metallicity is 



signifi cantly larger than in the case of H2BCG. This situation was first recognized bv lHenrv et al. 



(j2004l ) who dubbed it the sulfur anomaly, i.e., the systematic tendency for PNe to possess lower 
S abundances than do H II regions of similar oxygen abundance (metallicity) for reasons which 
have yet to be confirmed. Given our findings above concerning the similarities between PN and 
H2BCG abundance patterns for the alpha elements O, Ne, and Ar, we would expect S, another 
alpha element, to vary in lockstep with the others. 

An appreciation of the sulfur anomaly can be obtained by comparing the frequency distribu- 
tions of log(S/0), log(S/Ne), and log(S/Ar) for our PN sample and the H2BCG objects (Figs. [M 
18]) . The average of the logs along with the standard deviation is given in the legend in each case. 
In contrast to the large overlap visible in the comparisons of the other alpha elements above, there 
is a clear distinction between the two object groups in each of the three graphs. Comparing the 
averages of the logs of the abundance ratios we see that the offsets between the PNe and H2BCG 
are roughly 0.45 dex, with the PN average always falling substantially below that for H2BCG. In 
addition, the standard deviations indicate that the spread in each ratio is much greater for PNe 
than for H2BCG. 



In the original discussion of this situation in iHenry et al. I (120041 ) it was proposed that sulfur 



anomaly was likely related to problems inherent in the determination of the ICF for sulfur, since 
photoionization models indicate that in many PNe a significant amount of S is in the S"^'^ ionization 
stage, which is inaccessible at visual wavelengths. During the intervening time, however, numerous 
PNe have been studied spectroscopically in the IR by researchers using the Infrared Space Ob- 
servatory and Spitzer Space Telescope, both of which provide access to the [S IV] emission line 
centered at 10.5^m. These data obviate the need for using a sulfur ICF and provide for the direct 
computation of a total sulfur abundance simply by adding together the abundances of S+, S"*"^, 
and S+3. 



In Figs. [2] and [5] we have included the PN sulfur abundances (triangles) from lBernard-Salas et al 



(j2008l ) for a small sample of PNe whose S abundances were determined using IR data in the manner 
just described. One can readily see that the sulfur anomaly persists, as the scatter and system- 
atically low S abundances are similar to what we find in our own sample. Therefore, we can now 
state with some confidence that the sulfur anomaly is apparently unrelated to the sulfur ICF. 



It has been suggested by lPottasch and Bernard-Salas I (|2006l ) that the sulfur anomaly could be 



the result of dust formation. In particular, S may be removed by the formation of compounds such 
as MgS and FeS in those PNe exhibiting large S deficits. This wo uld be expected to occur more 



readi ly in C-rich environments, where sulfide formation is favored ([Lattimer et al. Ill978l : IZhang 



20091 ). Initial plots of the deviation of the S abundance from the expected value determined by 



the II2BCG track versus C/0 shows essentially no trend, where one might expect a correlation if 
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a larger C/0 is associated with C-rich environments where MgS can form. In addition, a simple 
calculation shows that the average S deviation would correspond to a drop of about 0.1 dex if in 
fact the sulfur deficit were due to MgS dust formation, an offset which is likely undetectable. 

5. Summary /Conclusions 

Continuing our previous work we present line strengths, abundances, and X/0 element ratios 
for a sample of 38 Galactic disk PNe bringing the total to > 120. The extended spectral coverage for 
this large sample of objects is atypical making available the [S III] lines at 9069A and 9532A for use 
in determining nebular diagnostics as well as total sulfur abundance. Including these objects in our 
compilation we've attempted to discern signatures of alpha element processing in PNe progenitor 
stars. We conclude the following: 

1. A careful comparison of the behavior of Ne, O, and Ar with respect to each other in PNe 
and II2BCG reveals systematically higher values of Ne/0 and lower values of Ar/Ne in PNe 
within the range of 8.5-9.0 of 12+log(0/II), suggesting that Ne is enhanced in many of these 
objects. The enhancement is likely coming from the extended processing of ^^N by alpha 
capture to yield ^^Ne. 

2. The positions of numerous outliers in several plots cannot be explained away by uncertainties 
in the abundance determinations. The progenitors of these objects may have formed from 
poorly mixed interstellar material or show significant synthesis or depletion of alpha elements. 

3. Type I and non-Type I PNe exhibit no difference in their alpha element abundance patterns. 

4. The sulfur anomaly, in which sulfur is found to be systematically lower in PNe than in H2BCG 
at the same metallcity, is reconfirmed. Recent work shows that the sulfur anomaly is likely 
not related to problems with the IGF and at first glance is too large to be accounted for by 
dust formation. 
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Table 1. Observation Information 



Object Offset {"y Total Blue Exp (sec) Total Red Exp (sec) Observatory^ 



Cn3-1 




429 


480 


K04 


H2-18 




1200 


1800 


K04 


Hb4 




1800 


1600 


K03 


Hb6 




800 


600 


K03 


He2-lll 




4260 


2940 


C04 


Hu 1-2 




1420 


1030 


K03, K04 


Ml-8 




4200 


3600 


C03 


Ml-13 




1620 


2300 


C03 


Ml-17 




960 


1200 


C03 


Ml-40 




1500 


1200 


K03 


Ml-42 




1650 


3360 


C04 


M2-52 




3600 


5000 


K03 


M2-55 




1800 


450 


K04 


M3-2 




3300 


3000 


C03 


M3-3 




3000 


3600 


C03 


M3-5 




2100 


2400 


C03 


Mel-1 




250 


1125 


K03 


Me2-2 




1400 


1080 


K03, K04 


Mz 2 




3180 


3612 


C04 


Mz 3 




1800 


2880 


C04 


NGC 1535 




2360 


2500 


C03 


NGC 2452 




1500 


1200 


C03 


NGC 5315 




961 


855 


C04 


NGC 6302 




700 


1350 


C04 


NGC 6369 


ION 


2000 


2200 


K03 


NGC 6445 


25N 


1200 


1200 


K03 


NGC 6537 




725 


300 


K03 


NGC 6741 




800 


400 


K03 


NGC 6751 


6S 


1500 


1500 


K03 


NGC 6778 




1500 


1800 


K03 


NGC 6803 




300 


150 


K03 


NGC 6804 


lOS 


1800 


5100 


K03 


NGC 6853 


83S, 120W 


2700 


1800 


K04 


NGC 6881 




1700 


500 


K03 


NGC 7008 


29N, HE 


1200 


1926 


K04 


NGC 7354 




3503 


4500 


K03 


PB6 




2000 


3000 


C03 


Vy 1-1 




360 


1440 


K04 



^From central star or center of nebula. 

= KPNO, C = CTIO, number denotes year of run. 



Table 2. Fluxes and Intensities 



Line 




Cn3-1 




H2-18 


Hb4 


Hb6 


He2-lll 


f(A) 


F(A) 


1(A) 


F(A) 


1(A) 


F(A) 


1(A) 


F(A) 


1(A) 


F(A) 


1(A) 


[O II] A3727 


0.292 


145 


165±37 


10.5 


24.9±5.55 


13.6 


37.9±8.47 


13.5 


48.7±10.93 


223 


416±95 


He II + Hll A3770 


0.280 


2.73 


3.10±0.68 


















He II + HIO A3797 


0.272 


3.36 


3.80±0.82 














2.87: 


5.15±1.85: 


He II + H9 A3835 


0.262 


5.15 


5.79±1.24 






2.61 


6.54±1.40 


1.70:: 


5.36±2.87:: 


4.25 


7.44±1.62 


[Ne III] A3869 


0.252 






41.2 


86.8±18.24 


38.4 


93.0±19.61 


31.5 


95.5±20.17 


150 


258±55 


He I + H8 A3889 


0.247 


10.8 


12.0±2.50 


4.36 


9.04±1.89 


5.40 


12.8±2.68 


3.90:: 


11.5±6.14:: 


16.2 


27.4±5.83 


[Ne III] A3968 


0.225 






10.4'» 


20.1±7.26=' 


12. 2=^ 


26.8±8.62^ 


9.96^ 


26.7±8.65=' 


42. 5=^ 


68.8±17.54'' 


He A3970 


0.224 


13.7 


15.1 


8.21=^ 


15. 9=^ 


7.27=^ 


16.0^ 


6.03^ 


16.2=^ 


10.4=^ 


16. 9=^ 


He I + He II A4026 


0.209 


0.557: 


0.612±0.211: 






0.887: 


1.85±U.d4: 






2.04:: 


3.19±1.69:: 


[S II] A4071 


0.196 


3.36: 


3.66±1.25: 






1.71: 


3.41±1.17: 


2.09:: 


4.95±2.61:: 


12.9 


19.6±3.87 


[S III] A9069 


-0.670 






99.7 


13.8±3.17 


510 


48.5±11.13 


471 


24.8±5.70 


256 


61.2±13.99 


P9 A9228 


-0.610 






34.6 


5.71±1.20 


37.2 


4.37±0.92 


49.2 


3.37±0.71 


6.74: 


1.83±0.64: 


[S III] A9532 


-0.632 


104 


78.4±17.08 


184 


28.5±6.19 


613 


66.5±14.45 


1941 


121±26 


849 


220±48 


P8 A9546 


-0.633 


8.03 


6.05±1.32 


12.0 


1.85±0.40 


13.1 


1.41±0.31 


61.0 


3.78±0.82 






c 




0.19 




1.28 


1.53 




1.91 




0.93 




Ha/H/3 




2.88 




2.i 


36 


2.87 




2.83 




2.78 




log Fh,?'' 




-10.90 




-12.96 




-12.10 




-12.22 




-12.88 





^Dcblcndcd. 

''ergs cm~^ in our extracted spectra 



Table 3. Temperatures and Densities 



Parameter^ 


Cn3-1 


H2-18 


Hb4 


Hb6 


He2-lll 


Valuc^ Notes'^ 


Value 


Notes 


Value Notes 


Value Notes 


Value Notes 


T[o III] 


8490±1134 


9601±438 




9241±408 


11180±1235 


15670±1162 


T[N II] 


7828±647 


10300 


Default. 


11720±1274 


12280±1167 


11770±876 


T[o II] 


7715±2605 


11610±5569 




15120±11090 


18160±14000 


9149±1969 


T[s II] 


7200±3504 






7612±3824 


11110±8984 


8828±2015 


T[s III] 


7206±435 Used 9532. 


10840±1808 


Used 9069. 


9988±874 Used 9069. 


11380±1075 Used 9532. 


15610±1978 Used 9532. 


Ne[s n] 


6195±5072 


4351±2840 




6001±4717 


4552±2952 


583±356 


Ne[ci ni] 


8451±5980 






5465±5439 


5367±1808 


1238±2726 



'^The emission lines (wavelengths in A) used to calculate electron temperatures and densities were: 5007, 4363 for T[o ni]! 6584, 5755 for T[n n]i; 3727, 
7324 for n]; 6716, 6731, 4071 for T[g „], 9532, 9069, 6312 for T[s m]; 6716, 6731 for Ne,[s n]; and 5518, 5538 for N^jci iii]- S 

'^Temperatures and densities given in kelvin and cm 
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'^If either 5755 or 6584 is unavailable T[n n] is estimated from T[o m] (|Kwitter fc Henry II2001I ). Given telluric effects when calculating T[g m] if the 
ratio of 9532/9069 is > 2.48 then 9532 is used as the nebular transition, if the 9532/90 69 ratio is less than 2.48 then 9069 is used. Default T[im uj and 
^c,[ci III] values and the high density limit for Nejs ii] are based on criteria discussed in iKwitter fc Henry I (|200l|). 



Table 4. Ionic Abundances 



Ion 




Cn3-1 




H2-18 




Hb4 




Hb6 




He2-lll 


T 


Abundance'' 


Tuscd 


Abundance 


Tused 


Abundance 


Tuscd 


Abundance 


Tuscd 


Abundance 


He+ 


[O III] 


5.57±0.70(-2) 


[O III] 


0.116±0.015 


[O III] 


0.121±0.015 


[O III] 


0.117±0.016 


[O III] 


0.144±0.021 


He+2 


[O III] 




[O III] 


3.82±1.23(-3) 


[O III] 


1.29±0.20(-2) 


[O III] 


1.97±0.30(-2) 


[O III] 


7.65±1.17(-2) 


icf(He) 




1.00 




1.00 




1.00 




1.00 




1.00 


O°(6300) 


[N II] 


*8.60±3.10(-6) 


[N II] 


*1.69±0.24(-6) 


[N II] 


*7.31±2.34(-6) 


[N II] 


*5.94±1.82(-6) 


[N II] 


*2.96±0.81(-5) 


00(6363) 


[N II] 


*8.22±2.98(-6) 


[N II] 


*8.24±4.21(-7) 


[N II] 


*8.23±2.65(-6) 


[N II] 


*5.87±2.46(-6) 


[N II] 


*2.93±0.81(-5) 


0° 


wm 


8.51±3.03(-6) 


wm 


1.57±0.23(-6) 


wm 


7.56±2.38(-6) 


wm 


5.92±1.83(-6) 


wm 


2.95±0.79(-5) 


0+{3727) 


[N II] 


*3.79±2.81{-4) 


[N II] 


•1.51±0.49(-5) 


[N II] 


*1.79±1.12(-5) 


[N II] 


*1.54±0.72(-5) 


[N II] 


*6.50±1.73(-5) 


0+{7325) 


[N II] 


*3.67±1.78(-4) 


[N II] 


*1.79±0.71(-5) 


[N II] 


*2.36±1.01(-5) 


[N II] 


*2.24±0.97(-5) 


[N II] 


*4.27±1.94(-5) 


0+ 


wm 


3.78±2.65(-4) 


wm 


1.65±0.33(-5) 


wm 


2.21±0.89(-5) 


wm 


2.12±0.85(-5) 


wm 


6.25±1.63(-5) 


0+2(5007) 


[O III] 


*1.35±0.71(-5) 


[O III] 


•5.35±1.27(-4) 


[O III] 


*5.78±1.36(-4) 


[O III] 


*3.67±1.46(-4) 


[O III] 


n.l8±0.33(-4) 


0+2(4959) 


[O III] 


*1.32±0.67(-5) 


[O III] 


*5.06±0.96(-4) 


[O III] 


*5.51±1.04(-4) 


[O III] 


*3.54±1.30(-4) 


[O III] 


*1.16±0.26(-4) 


0+2(4363) 


[O III] 


*1.35±0.71(-5) 


[O III] 


*5.35±1.27(-4) 


[O III] 


*5.78±1.36(-4) 


[O III] 


*3.67±1.46(-4) 


[O III] 


n.l8±0.33(-4) 


0+2 


wm 


1.34±0.70(-5) 


wm 


5.28±1.16(-4) 


wm 


5.71±1.25(-4) 


wm 


3.64±1.41(-4) 


wm 


1.18±0.30(-4) 


icf(O) 




1.00 




1.03 




1.11 




1.17 




1.53 


S+ 


[N II] 


*1.47±0.92(-6) 


[N II] 


*1.03±0.34(-7) 


[N II] 


*5.82±3.23(-7) 


[N II] 


*4.71±2.03(-7) 


[N II] 


*2.80±0.67(-6) 


S+(6716) 


[N II] 


1.48±0.92(-6) 


[N II] 


1.03±0.34(-7) 


[N II] 


5.80±3.26(-7) 


[N II] 


4.70±2.03(-7) 


[N II] 


2.83±0.68(-6) 


S+(6731) 


[N II] 


1.47±0.92(-6) 


[N II] 


1.03±0.34(-7) 


[N II] 


5.84±3.22(-7) 


[N II] 


4.72±2.02(-7) 


[N II] 


2.78±0.66(-6) 


S+ 


[S II] 


1.80±3.01(-6) 


[S II] 




[S II] 


1.28±2.04(-6) 


[S II] 


5.40±7.72(-7) 


[S II] 


4.56±2.77(-6) 


S+2(9532) 


[S III] 


*7.55±2.51{-6) 


[S III] 


•1.77±0.62(-6) 


[S III] 


*6.96±2.27(-6) 


[S III] 


*5.23±1.72(-6) 


[S III] 


•5.20±1.67(-6) 


S+2(6312) 


[S III] 


*7.55±2.51(-6) 


[S III] 


*1.77±0.62(-6) 


[S III] 


*6.96±2.27(-6) 


[S III] 


*5.23±1.72(-6) 


[S III] 


*5.20±1.66(-6) 




wm 


7.55±2.51(-6) 


wm 


1.77±0.62(-6) 


wm 


6.96±2.27(-6) 


wm 


5.23±1.72(-6) 


wm 


5.20±1.67(-6) 


icf(S) 




1.00 




1.83 




1.81 




1.68 




1.17 



''wm indicates the calculated mean value for the ionic abundance weighted by the observed flux 
''an * indicates those values included in the weighted mean 



Table 5. Total Elemental Abundances 



Parameter Cn3-1 H2-18 Hb4 Hb6 He2-lll Solar Ref ^ Orion Ref 



rie/ rl 


0.0 ( =tu. 






U. izUztU, 


.UiO 




n 1 Q/i_i_n 


.UiO 




U.ioDztU, 


.Ul ( 




U.zziztu 


.UZD 




o.oii^-z ) 


y 


.oyJ[-z ) 


N/H 


7.75±2, 


.50(- 


■5) 


6.24±1, 


.93(- 


■5) 


5.17±1, 


.36(- 


-4) 


4.09±1, 


.66(- 


■4) 


8.25±2, 


.23(- 


-4) 


6.03(-5) 


6 


.03(-5) 


N/0 


0.198±0 


.096 




O.llliO, 


.025 




0.787±0, 


.133 




0.907±0, 


.181 




2.99±0, 


.70 




0.132 




0.115 


0/H 


3.91±2 


.66(- 


-4) 


5.62±1, 


.20(- 


-4) 


6.57±1, 


.40(- 


-4) 


4.50±1, 


.65(- 


■4) 


2.76±0, 


.62(- 


-4) 


4.57(-4) 


5, 


.25(-4) 


Ne/H 








1.12±0, 


.26(- 


-4) 


1.51±0, 


.35(- 


-4) 


7.76±3 


.25(- 


■5) 


1.27±0, 


.34(- 


-4) 


6.92(-5) 


7, 


.76(-5) 


Nc/0 








0.199±0, 


.032 




0.230±0 


.038 




0.172±0, 


.030 




0.461±0 


.079 




0.151 




0.148 


S/H 


9.02±2, 


.89(- 


-6) 


3.41±1, 


.39(- 


■6) 


1.37±0, 


.52(- 


■5) 


9.60±3, 


.70(- 


■6) 


9.36±2, 


.61(- 


-6) 


1.38(-5) 


1, 


.48(-5) 


S/0 


2.31±1, 


.51(- 


-2) 


6.07±2, 


.43(- 


■3) 


2.08±0, 


.81(- 


■2) 


2.13±0, 


.78(- 


■2) 


3.39±1 


.05(- 


■2) 


3.02(-2) 


2 


.82(-2) 


Cl/H 


2.16±0, 


.57(- 


■7) 








1.30±0, 


.41(- 


■7) 


1.07±0, 


.25(- 


■7) 


1.93±0 


.58(- 


■7) 


3.16(-7) 


2, 


.14(-7) 


Cl/0 


5.53±3, 


.40(- 


-4) 








1.97±0, 


.67(- 


-4) 


2.38±0, 


.76(- 


-4) 


7.00±2, 


.22(- 


-4) 


6.92(-4) 


4, 


.07(-4) 


Ar/H 


6.76±2, 


.68(- 


■7) 


1.46±0, 


.29(- 


■6) 


3.91±0 


■71(- 


-6) 


4.05±1, 


.14(- 


■6) 


4.75±0, 


.85(- 


-6) 


1.51(-6) 


3 


.09(-6) 


Ar/0 


1.73±L 


.29(- 


■3) 


2.60±0, 


.48(- 


■3) 


5.96±1, 


.04(- 


■3) 


9.00±1, 


.56(- 


■3) 


1.72±0, 


.28(- 


■2) 


3.31(-3) 


5, 


.89(-3) 



''all abundances are of the form N(X) 

''Solar values taken from Asplund et al. 2005. 

'^Orion values taken from Esteban et al. 1998 



-22- 



Table 6. Linear Fit Parameters 



Relation^ Y Intercept Slope Correlation Coefficient Objects 



H II Regions & BCGs 



Ne/H V. 0/H 


-0.84±.27 


+1.02±.03 


+0.96 


85 


S/H V. 0/H 


-1.63±.22 


+1.01±.03 


+0.97 


73 


Ar/H V. 0/H 


-3.03±.75 


+1.10±.09 


+0.80 


82 


S/H V. Ne/H 


-0.46±.26 


+0.95±.04 


+0.96 


69(2) 


Ar/H V. Ne/H 


-1.26±.30 


+0.96±.04 


+0.94 


74 


S/H V. Ar/H 


+0.90±.19 


+0.97±.03 


+0.96 


69 


Ne/0 V. 0/H 


-0.84±.27 


+0.017±.03 


+0.055 


85 


Ar/Ne v. Ne/H 


-1.26±.30 


-0.036±.04 


-0.10 


74 


Planetary Nebulae 


Ne/H V. 0/H 


-1.06±.55 


+1.05±.06 


+0.81 


146 


Cl/H V. 0/H 


-3.66±.73 


+1.01±.08 


+0.72 


134 


Ar/H V. 0/H 


-1.76±.68 


+0.94±.08 


+0.70 


153 


Cl/H V. Ne/H 


+0.38±.71 


+0.59±.09 


+0.51 


128 


Ar/H V. Ne/H 


-1.77±.47 


+0.57±.06 


+0.63 


147 


Cl/H V. Ar/H 


-0.35±.40 


+0.85±.06 


+0.77 


132 



^X/H values are expressed as 12+log(X/H), while X/Y values are given as 
log(X/Y), where X and Y are metals. 



Table 7. Model Results: Yields from Massive Stars and SNIa 



Relation Predicted Slope Observed Slope 



Ne/H V. 0/H +1.11 +1.02 

S/H V. 0/H +1.20 +1.01 

Ar/H V. 0/H +1.06 +1.10 

S/H V. Ne/H +1.10 +0.95 

Ar/H V. Ne/H +0.96 +0.96 

S/H V. Ar/H +1.13 +0.97 
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Fig. 1.— 12+log(Ne/H) versus 12+log(0/H) for PN data (HKB04, Milmgo09, and Kwitter09) 
along with H II region and blue compact galaxy abundance data (H2BCG) and the Sun. References 
are provided in the text. The solid line is a least squares fit to the H2BCG data. 
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Fig. 2. — Same as Fig. [T]but for S versus O. 
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Fig. 3. — Same as Fig. [T]but for CI versus O. No CI data are available for the H2BCG objects. 
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Fig. 5. — Same as Fig. [T]but for S versus Ne. 
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Fig. 6. — Same as Fig. [3] but for CI versus Ne. 
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Fig. 7. — Same as Fig. [T] but for Ar versus Ne. 
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Fig. 8. — Same as Fig. [3] but for CI versus S. 
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Fig. 9. — Same as Fig. [T] but for S versus Ar. 
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Fig. 10. — Same as Fig. [3] but for CI versus Ar. 
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— MWG Disk: -0.64 (.138) 
■ ■ LD:SMC -0.86 (.25) 

- ■ H2BCG: -0.70(.11) 
-- LD:LMC -0.77 (.25) 




log(Ne/0) 



Fig. 11. — Frequency di stribution of log(Ne/0) fo r the Galactic PNe in this study, the PNe in 
the LMC and SMC from lLeisv fc Dennefeldl (|20od ). and for H2BCGs. The numbers in the legend 
represent the average of the logs (and the standard deviation). 



20 



5 





1 ' 1 




' 1 ' 


1 ' 1 






— PNe: -1.66(.27) 

- H2BCG:-1.52(.11) 


, I 
i |i 

/ 'i 












/ ■' 
/ ■' ' 

!?' i 








1 /^a/i 




1 


i\a a 1 





log(Ar/Ne) 

Fig. 12.— Same as Fig. [U] but for log(Ar/Ne) 





Fig. 14. — Frequency distribution of log(Ne/0) differential offsets for Type I PNe and non-Type 
I PNe in our sample. The numbers in the legend indicate the average of the log values and the 
standard deviations. 
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Fig. 16. — Frequency distribution of log(S/0) for our PN sample and H2BCG. Averages of the logs 
and standard deviations appear in the legend. 
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Fig. 17. — Same as Fig. [TBI but for log(S/Ne). 




Fig. 18.— 



Same as Fig. [TBI but for log(S/Ar). 



